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Abstract. In this paper, we propose a strategy to determine the Dirichlet-to-Neumann (DtN) 
operator for infinite, lossy and locally perturbed hexagonal periodic media. We obtain a factoriza- 
tion of this operator involving two non local operators. The first one is a DtN type operator and 
corresponds to a half-space problem. The second one is a Dirichlet-to-Dirichlet (DtD) type operator 
related to the symmetry properties of the problem. The half-space DtN operator is characterized via 
Floquet-Bloch transform, a family of elementary strip problems and a family of stationary Riccati 
equations. The DtD operator is the solution of an affine operator valued equation which can be 
reformulated as a non standard integral equation. 
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1. Introduction and problem setting 

Periodic media appear in many physical and engineering applications related to 
wave type systems. In solid state theory, it is well-known that the existence of in- 
sulators and conducting materials can be explained by the periodic properties of the 
crystal. More precisely, within the one-electron model of solids, electrons move under 
a Hamiltonian with periodic potential yielding the existence of band gaps (see for ex- 
ample Reed and Simon |34l p. 312]). In optics, many devices used in microtechnology 
and nanotcchnology involve materials with such electromagnetic properties, known as 
photonic crystals. For a general introduction to the physics of photonic crystals, we 
refer the interested reader to the monographs by Joannopoulos et al. [23] , Johnson 
and Joannopoulos Sakoda [37] or the review paper by Busch et al. 8]. Con- 
cerning mathematical aspects related to photonic crystals, see for instance the book 
of Kuchment |28| . About twenty years ago, the analysis of elastic wave propagation 
in periodically structured media lead to the concept of phononic crystals (see for in- 
stance the book by Maldovan and Thomas [3"T]). 

Although they concern materials involving quite different scales and types of waves 
(acoustic, electromagnetic, elastic), photonic, phononic or real crystals (in solids) 
share a lot. The main common feature, inherited from the periodic structure, is the 
appearance of band-gaps (forbidden bands), i.e. the strong attenuation of a certain 
range of frequencies (at least in some directions). This phenomenon is due to the fact 
that an incident wave on the crystal is multiply scattered by the periodic structure, 
leading to possibly destructive interference (depending on the characteristics of the 
crystal and the wave frequency). 
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Most applications involving periodic structures use these band gaps to control wave 
propagation properties of materials (see [5] for an extended list of references). More 
precisely, the aim can be to prohibit propagation of waves [35], allow propagation of 
some frequencies and/or some directions like in optical filters [33], localize waves by 
creating a defect in the periodic structure [301 EZ] , create optical nanocavities p] or 
to guide waves with bends [35], Y-j unctions [B] and T-junctions [T4"] . 
In the above applications, numerical simulation plays a crucial role in the design of 
photonic, phononic and real crystals. In particular, the problem of determining which 
physical and geometrical properties of a perfect or imperfect crystal yield band gaps 
at specified frequencies, guiding waves along some desired paths or trapping them in a 
given cavity is computationally challenging. The complexity of designing such devices 
is increased by the strong scattering effects involved, due to the high contrast of the 
composing materials on the one hand, and the geometric scales (wavelength and size 
of the periodic cell) on the other hand. Therefore, several numerical methods and 
tools have been developed to simulate wave propagation in infinite periodic media. 
A first class of methods covers problems where the periodicity can be handled via 
homogenization techniques [5J [7] . Indeed, when the wavelength is much larger than 
the period, the periodic medium behaves asymptotically like a homogeneous one. The 
unboundedness of the medium can then be treated analytically using integral equa- 
tions, Dirichlet-to-Neumann (DtN) operators or Perfectly Matched Layers techniques. 
On the contrary, a second class of methods keeps the periodicity but considers only 

• media which are finite (i.e. constituted of a finite number of periodic cells 
embedded in an infinite homogeneous domains) [HI [I3J [53] [33] ; 

• media which can be reduced to finite domains. The case of imperfeclQperiodic 
media has been more intensively investigated [13 E3 HI] ■ In these works, the 
problems considered lead to exponentially decaying solutions, which allows 
the truncation of the infinite periodic media. 

More recently, DtN operators have been derived for infinite two-dimensional peri- 
odic media containing local defects in [171 119] . The main assumptions concern the 
directions of periodicity (orthogonal), the corresponding periodicity lengths (commen- 
surate) and the presence of a dissipative term (arbitrarily small). 



In this work, we are interested in the analysis of hexagonal periodic media containing 
a local defect (see Figure [TTTj ) . We call hexagonal periodic medium (also known as 
hexagonal lattice) a two dimensional domain where 

• the angle between the periodicity directions is 7r/3; 

• the periodicity cell has hexagonal symmetry (see Definition 3.1 1. 

Note that like in (TTl [19] , one could consider the media as a periodic one with two 
orthogonal directions of periodicity. However, the corresponding periods would not 
be commensurate in this case. 

Hexagonal lattices appear in quantum mechanics [10l E3 [34] , phononics [29] 321 EE] 
and photonics [HI EJ [20] . Although they concern different applications and involve 
different types of waves, the corresponding problems are quite similar from the math- 
ematical point of view. Essentially, as far as numerical simulation is concerned, the 
main issue is to determine a transparent boundary condition to reduce the problem 
initially set on an unbounded domain to a problem set on a bounded one contain- 



1 A local defect is generally created by changing locally the physical/geometrical properties of the 
crystal. 



Transparent boundary conditions for locally perturbed infinite hexagonal periodic media 3 




oooooooo 
ooooooooo 
oooooooo 



ooooooooo 

AVAV ViVwJ 0000 

Mi 000 



OOOOO0O0 

OOOOOO <3><3 O 

oooooooo 



Fig. 1.1. A locally perturbed photonic crystal with triangular lattice (cross section of a photonic 
crystal optical fiber (left) and representation of the corresponding hexagonal periodicity cell (right)). 



ing the defect cell. The differential operators involved are the ones describing the 
underlying physics of the problem under consideration, namely 

• Quantum mechanics: in a (classical) crystal, the mathematical formulation 
of the problem leads to the Schrodinger operator 

Au := -Aw + (V + ip)u (1.1) 



where V(x) denotes the potential in the lattice and p is the Laplace variable. 
Phononics: the operator involved is the elasticity system (see [31 HI 1251 

gehsehui) 

Au := -diver(u) + w 2 pu, (1.2) 

where u denotes the displacement and er(u) — C : e(u) the stress tensor, in 
which e(u) stands for the strain tensor and C the 4t/i-order elasticity tensor 
(/? and lo respectively denote the mass density and the time frequency). 
Photonics: in this case, electromagnetic propagation is described by the 
vector Maxwell's equations. For two-dimensional photonic crystals, these 
equations reduce to the following scalar equations respectively in the cases of 
TE (Transverse Electric) and TM (Transverse Magnetic) polarizations: 

TE case: Au := Au + uj 2 n 2 u 
TM case: Au := -div (-^ Vit^ + uj 2 u 

where n(x) denotes the index of refraction and ui the wavenumber. 
Note that all the above operators involve an elliptic principal part. Moreover, the 
time dependence has been eliminated via a Fourier (or Laplace) transform, leading 
to the appearance of the real parameter lu (or the complex parameter p). Perfect 



periodic media are described via the operators (1.1), (1.2) or (1.3) -depending on 
the considered application- with periodic coefficients (potential, elasticity tensor and 
mass density or index, respectively for classical, phononic and photonic crystals). 



The introduction of a defect in the periodic structure is taken into account by adding 
a bounded obstacle or locally perturbing the coefficients. Typically, scattering by an 
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impurity in a classical crystal (i.e. when the lattice contains a different atom at one 
of the lattice points), we have 

V = Vpcr + V 

where V pcl is periodic and Vo is a short range potential describing the local perturba- 
tion (see 01 p. 312]). 



df 




Fig. 1.2. The hexagonal periodic medium with defect. 



In this paper, we will restrict our study, for the sake of simplicity, to the case of an 
infinite photonic crystal f2 = R 2 containing a localized defect in the hexagonal cell fl l 
(whose side length is denoted by d). We denote by E l = the boundary of this 
cell and by il e = \ il l its exterior. From the mathematical point of view, the model 
problem we consider throughout the paper is a dissipative Hclmholtz equation 

Au + pu = f, infi, (1.4) 

where the following assumptions will always be supposed to hold true 

• (Al) p is a local perturbation of a hexagonal periodic function p pcr . More 
precisely: 



Ppc 



- Po 



where 



- for all x = (x, y) <E and all (p, q) G Z 2 , p peI (x + pei + qe 2 ) = p pe r(x) 
with ei = (Zd/2, V3d/2) and e2 = (0, yf&d) the two directions of peri- 
odicity of the media (see Figure 1.2); 



— Pper and po have hexagonal symmetry (see Definition 3.2); 

- Supp(p ) C fP. 

(A2) p satisfies the dissipation property 

|Imp(x)| > p b > 0, VxeO. (1.5) 
(A3) The source / is compactly supported in SI* and has hexagonal symmetry. 
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Before going further, let us make some comments about the above assumptions. In 
(A2), condition (1.5| guarantees the existence and uniqueness of finite energy solu- 
tions of (1.4) (i.e. solutions in in the infinite domain SI (as it can be easily 
checked by using Lax-Milgram Lemma). When no dissipation is assumed (i.e. when 
p is real- valued) , such existence and uniqueness issues will not be discussed in this 
paper, nor will be the derivation of a limiting absorption principle. Let us emphasize 
that the latter remains, to our knowledge, an open question (see |19l Remark 2] and 
[2"5] for a similar discussion). 



The main goal of this paper is to propose a method to solve (1.4) in the infinite 
domain fl under assumptions (A1)-(A2)-(A3). The main steps of our approach are 
presented in a formal way in the next Section. 

2. Formal presentation of the method For the reader's convenience, we 
describe in this section the broad outlines of our approach in order to help the reader 
getting an overview of the proposed strategy. Willfully, we decided to focus on the 
description of the mains steps of the method and to skip in this part all the technical 
details related to the functional framework. 

Our approach is adapted from the one used in |17j for the case of square lattices. 



The key idea is to reduce problem (1.4), which is set in the unbounded domain f2, to 
a boundary value problem set in the cell SI 1 containing the defect. To achieve this, 
we need to derive a suitable transparent boundary condition on E 1 associated with a 
DtN operator A. More precisely, we note that the restriction u l :— u|fj» & H 1 ^ 1 ) 
solves the interior boundary value problem 

Au { +/»«* = /, mil 1 , 
du 1 

^-+Au* = 0, on£\ 
where v % is the outgoing unit normal to Jl l and A denotes the DtN operator defined 

by 



(2.1) 



in which u e (<fi) G H (Q e ) is the unique solution of the exterior problem 

f Au e (cf>)+pu e (4>) =0, inO e , 
\ u e {cf>) = <f>, on £\ 

The core of the paper is thus devoted to the computation of this DtN operator A. 
First of all, let us emphasize that due to the symmetry properties of the original 
problem (27r/3 rotational invariance, see Assumptions (Al) and (A3)), it suffices to 
compute this DtN operator A for Dirichlet data cf> on S 1 with hexagonal symmetry 
(See section [3] for more details). In this Section, we restrict our analysis from now on 
to such symmetric data. 

In order to describe the different steps of our method, we need to introduce some 



additional notation. Let Yj H be the boundary depicted in Figure 2.1 and let il H 
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Fig. 2.1. The half-space DtN operator A 1 





Fig. 2.2. The DtD operator D 2n / 3 



be the half-space to right of Y, H . Given a Dirichlet data <f> on ^ H \ we define the 
half-space DtN operator by setting (see Figure 2.1 1 



A 



(2.3) 



where v H is the outgoing unit normal to Sl H and u H ((f>) is the unique solution in 
ff 1 (A,il ff ) of the half-space problem 



Au H ((f)) + pu H ((/)) =0, 
u H {<p) = <f>, 



in 

on E H . 



(2.4) 



Transparent boundary conditions for locally perturbed infinite hexagonal periodic media 7 



The second important ingredie nt ne eded is the DtD operator Z? 27r /3 defined from the 
boundary Yi 1 to E ff (see Figure 2.2) by the formula: 



(2.5) 



where u e (<j)) is the unique solution of the exterior problem (2.2) for a given Dirichlet 
data on S'. 

Step 1 : Factorization of the DtN operator A. With the above notation, it 
is clear that given a Dirichlet data <f> on T, 1 , the functions u e ((j>) and u H (D 27r /y,4>) are 
both solutions of the half-space homogeneous Helmholtz problem 

AU + pU = 0, 

with the same Dirichlet condition on namely D 2n / 3 cf>. The uniqueness of the 
solution of this problem implies that 

u e ((f>) |n* = u H (D 2 „/ 3 <j>), 

and in particular, the corresponding normal derivatives coincide on the part of the 
boundary D Y> H where they are both defined, yielding 



(2.6) 



the remaining part of A</> on S l being recovered using hexagonal symmetry. The 
above relation constitutes the starting point of our strategy : the DtN operator A 

a 



can be computed via the factorization formula (2.6). We provide in Theorem 



4.1 



more precise statement of this factorization result, paying a particular attention to 
the functional framework. Thus, our problem is now reduced to the computation of 
the half-space DtN operator A H and the DtD operator D 27r /3 for well prepared data 
(i.e. having 2tt/3 rotational invariance). 

Step 2 : Characterization of the DtN operator A H . In order to compute 
the half-space DtN operator A H , the first key ingredient is the (partial) Floquet- 



Bloch transform in the vertical direction (see § 5.1 1 . More precisely, applying it to the 



Helmholtz half-space problem (2.4), we will see that it suffices to consider the case of 



k— quasiperiodic Dirichlet data (f> on that is 

0(y + qL) = 0( y y< lkL , VyeM, VgeZ, 

where L = \/3d denotes the period in the vertical direction and k € (—tt/L,tt/L). 
From now on we will thus restrict our analysis to k— quasiperiodic Dirichlet data. 



The second im portan t tool we need is the so-called propagation operator Vk, defined 
as follows (see { 5.2.1 ). For any k— quasiperiodic Dirichlet data <j> on Y, H , Vk 4> is noth- 
ing but the trace of the solution u H (<f>) of (2.4) on the trasnalted vertical boundary 
Y, H = Yi H + ei (see Assumption (Al)). The main advantage of this operator is that 
it allows us to determine the solution u ((f)) of the half-space problem in any cell from 
the knowledge of u H ((f>) on a reference cell. More precisely, with the notation given 



in Figure 2.3 we can prove thanks to a uniqueness argument that 



u H (4>)\ 



JqkL H 



u H ((V k )^ 
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d 




Fig. 2.3. Description of the half-space Q and its periodicity cells 



According to the above relation, solving the half-space problem amounts to determin- 
ing the propagator Vk and the values of u H (</>) but only in the reference cell Coo- F° r 
the latter, by linearity, this can be achieved by solving two elementary cell problems 
set in Coo (see § 5.2.2). Regarding the determination of the propagator operator Vk, 

Writing the matching of the normal derivatives of 



5.2.3) 



we proceed as follows (see 
the solution u H ((f>) across a suitable part of the interface we will show that Vk 
solves a stationary Riccati equation. 

Step 3 : Characterization of DtD operator D 27T /3- While the determination 
of the half-space DtN operator A H described in Step 2 essentially uses the periodicity 
properties of the medium, the determination of the DtD operator -D 27r /3 crucially uses 
its hexagonal symmetry. More precisely, we show in Theorem 



6.1 



that D 2 n/3 solves an 



affine operator- valued equation (see equation (6.1 1 in Theorem 6T I that is well-posed. 
In order to handle this equation from a practical point of view, we use once again 
Floquet-Bloch variables instead of the physical ones. Doing so, we reduce this affine 
equation to the solution of a set of non standard constrained integral equations (see 



equations (6.9)). 

Remark 2.1. For the sake of clarity, we have restricted our analysis in this paper 
to the model problem (1.4), under assumptions (A1)-(A2)-(A3). From the physical 
point of view, this problem describes the radiation in a photonic crystal of a source 
localized ( typically a real source ) in one cell of the periodic medium. We would like 
to emphasize that our approach can be generalized to tackle other types of problems: 
extended defects (i.e. covering more than one periodic cell), non penetrable defects or 
scattering problems by a local inhomogeneity / 'obstacle in the crystal. Non symmetric 
configurations can also be considered but involves more intrincate integral equations 
(see \1 7| / for non symmetric square lattice). 

The precise description of the three steps briefly outlined above involves many 
technical details related to the functional framework that needs to be carefully set- 
tled for the DtN and DtD operators. In Section [3j we collect some notational and 
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mathematical background related to geometries with hexagonal symmetry. Section [4] 



provides the factorization result (Theorem 4.1 1, which plays a key role in our method 



for computing the DtN operator A. The next two sections are devoted to the analysis 
of the two non local operators A H and -D 27r /3 involved in the factorization. More pre- 
cisely, Section [H] deals with the half-space DtN operator A H , while the DtD operator 
^2tt/3 is studied in Section [6j An algorithm will be given in conclusion to summarize 
this construction. 

3. Hexagonal symmetry and related results 

3.1. Definitions and functional framework 

Definition 3.1. A domain O ofM. 2 has hexagonal symmetry if there exists a rotation 
of angle 2tt/3, denoted 62^/3 for which O is invariant. In the above definition, the 
center of the rotation is not specified to simplify the notation. In the following, this 
center will always be the center of gravity of O. 

Remark 3.1. If a domain O ofM 2 has hexagonal symmetry, its boundary dO also 
has hexagonal symmetry. 

Definition 3.2. Let O be an open set with hexagonal symmetry and let g be a real 
or complex valued function defined on O. Then, g has hexagonal symmetry if 

9=9° @2tt/3- 

Let O be an open domain with hexagonal symmetry (typically O = Q e , f2 J and 
dO = For s > 0, we denote by H%. 3 (0) the closed subspace of H'(0) defined 

by 

H^ /3 (0) = {ve H s (0), v = vo 9 27r/3 }. (3.1) 
One can easily check that the restriction to H^^iO) of the trace operator 70 defined 

by 

Vu G H\0), 1q u = u\ do g H 1/2 {dO), 
defines a continuous operator from H 27T ^ 3 (0) onto H^ 3 (dO). 

Finally, let us define the appropriate functional space for the normal trace of a function 
with hexagonal symmetry. To achieve this, we introduce the spaces 

H x (A,0) = {ue H\0) I Au G L 2 (0)} 

and 

Hl /3 (A,0) = {ue Hl /3 (0) I Au G L 2 2rr/3 (0)}. 

The next result expresses in term of functional spaces the commutativity of the 
Laplace operator with the rotation 62^/3 ■ 

Theorem 3.3. Let O be an open set ofM 2 with hexagonal symmetry. The Laplace 
operator commutes with any unitary transform and maps H^i 3 (A, O) onto L\ v / 3 (0). 

We can now extend Definition (3.1 ) to the space H~ 1 ^ 2 {dO) where O is with hexagonal 



symmetry (typically O = fl e or Q l and dO = E 4 ) 

Definition 3.4. Let O be an open set of M 2 with hexagonal symmetry. We define 
the closed subspace H~\^(dO) of H-^ 2 {dO) by 



du 

dv 



,ueHl /3 (A,0)} 
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where v is the outgoing unit normal to O. 

Obviously, 71 is a continuous application from H 27T ^ 3 (A,0) onto H 2 ^j^{dO). 

3.2. Restriction and extension operators Let R be the restriction operator 
denned by 

R : L 2 (E J ) -> L 2 (E°) 



where E° is the right part of E* (see Figure 1.2). One can easily check that R defines 
an isomorphism from the subspace L\^, S (Y?) := {v G L 2 (E l ), v — v o 62^/3} onto 
L 2 (E°) that we shall denote by R 27T /3- Its inverse E 27T / 3 is an extension operator 
which can be given explicitly thanks to the rotation 6)2^/3 : 





E 27r / 3 4> 






6L 2 (S°), 


E^/30 


e 27r/3 £° = ' 


i>° ©-27T/3 






©L /3 E ° = ' 





We define the following space 

H li%(z°) ■= {R2,/ 3 <P,^eH 1 2 i%(^)} ) 

= {^^(E ), i^G^^E*)}. 

This space is nothing but the space of functions in iJ 1 / 2 (S°) which are periodic. 

^ t 2 

We now explain how to extend the restriction operator R 2jv /3 to H 2 ^L (E 1 ) (see 
Definition 3.4). This can be done by duality, noticing that we have 

g £L/ 3 (^), G L 2 (£°), (i2 27r/3 0,V) EO = I (0,^/3^)^ , 
where (•, -)s° (resp. (v)e*) is the scalar product in L 2 (E°) (resp. L 2 (Y> 1 )). This last 

— 1/2 

relation suggests an extension of R 27T /3 to H 2 ^L (£ l ) by 



v <^ e H £/s(F), G ^ 3 (E U ), (i^/a^EO - ^^W/s^s*, 

where (•, -)so (resp. (•, •)£*) is the duality product between the two spaces [H 2 ^ 2 , 3 (TP)]' 
and H 2 ^ 2 , 3 (S°) (resp. H~^j 3 2 (Y, 1 ) and iJ^^E 1 ) ). We introduce the closed subspace 



1-1/2 , v0 



of 



H. 



-l/2/ v 



2tt/3 



and conclude that R^/3 is a linear continuous map from i? 2 J/g (£*) onto H 2 ^j 2 (E°). 

Analogously, E 27T /3 can be extended to a linear continuous mapping from ff 2 V 3 2 (E°) 
onto H~^(T. i ) using: 

V ^ G H^p ), V0 G ff^E*), (^/3^ 0>E« = 3(V, i?2./30) S o. 
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3.3. Symmetry properties for the exterior problem Using the definitions 
and properties of media and functions with hexagonal symmetries and considering 
assumptions (Al) and (A2), one can show the following theorem: 

1/2 \ 1 

Theorem 3.5. If cf> € iJ 2 ^y 3 (S J ) ; then the unique solution u e (4>) of (2.2) belongs to 

^271-/3(^1 an d A</> € ^2n/3 > w ^ ere A is the DtN operator defined in (2.1 1. 

Proof. It suffices to prove that u e (0 2lr /30) is also a solution of (2.2 1 and then 
conclude using a uniqueness argument. The last implication follows immediately from 
the definition of the norm al tr ace operator. □ 

A consequence of Theorem 3.5 is that A maps continuously H^^{Yj 1 ) onto H 2 tt/3 (^ l )> 
leading to the natural definition of the DtN operator 



A 



2tt/3 



1/2, 



2tt/3 



t/ 



/3 



(3.2) 



Finally, taking into account the assumptions (A1)-(A3) concerning the hexagonal 
symmetry of the local perturbation po and of the source term /, and using uniqueness 
argument, we deduce that 

Theorem 3.6. Let u be the unique solution of (1.4). Then the restriction u 1 — u\ Qi 
of u is the unique solution of the interior boundary value problem 

A^+ptf = /, intt 1 , 

— - + A 2v/3 u l = 0, onE 1 . 
ov l 1 

The rest of the paper is devoted to the determination of the DtN operator A 27r /3. 

4. Factorization of the DtN operator 

First o f all let us recall some useful notation. Let Yi H be the boundary depicted in 
Figure 2.1 and let fl H be the half-space to right of T, H . 



Let A H £ CiH 1 / 2 ^")^- 1 / 2 ^)) be the half-space DtN operator (see Figure |2.1) 

A**=*£<£ (4.1) 



H 



where v H is the exterior normal to ft H and u H (4>) is the unique solution in if 1 (A, H, H ) 
of the half-space problem 



Au H ((j)) + pu H ((l))=0, 
u H (4>) = 0, 



in n H , 

on E H . 



(4.2) 



Let L> 27r / 3 e CiH 1 / 2 ^),!! 1 / 2 ^ 11 )) be the DtD operator defined by (see Figure 2.2): 



(4.3) 



where u e (<f>) is the unique solution of (2.2). 



Moreover, we need to introduce a restriction operator R H from H H to S°. As we 
need to apply this operator to functions of _ff~ 1 / 2 (E ff ), R H has to be defined in a 
weak sense. We denote by E H G £(L 2 (S°), L 2 (T, H )) the extension operator by zero 
from E° to Y} H : 



E H d> = 







on S°, 

on T, H \ E°. 
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Let #o /2 (£ ) be the subspace of H^ 2 (T, ) defined by: 

i4 /2 (£°) = {0 e 1/2 (E°) | E H <t> G ff 1/2 (E ff )}, 

and let £r _1 / 2 (E ) = ^^ 2 (E )^ be its dual space. The restriction operator R H G 
£(i7- 1 /2( S H) ; ^-i/2( E 0)) can thcn be defined by duality 

(R H (j>, V , ) i g--l/2( E 0), fl -l/ o 2 ( S 0) = (4>1 E H lp) H -l/2( T ;Hy H l/2( S H^ 

V(0, V) G H-Va( E « ) x ff# 2 (£°). (4.4) 
The main result of this section reads as follows. 

" > "2tt/3 



Theorem 4.1. TTie operator R H o A H o D 2lT / 3 maps H^ 2 /3 (T, i ) into H 2 Y 3 2 (£ ) and 



A 27r /3 G £( J EZ^^ 3 (E*),iZ 2lr 1 / / ^ 2 (E*)) defined by (3.2 I admits the factorization 



e HV 2 /3 (p), A 27r/3 = £ 27r/3 oij fl oA fl o £> 2V3 </> 



where 



D 2n/3 G £( J ff 1 / 2 (E i ),ff 1 / 2 (E ff )) is the DtD operator defined by Q 



A ff G £(i7 1 / 2 (E H ),i/- 1 /2( S ^)) i s i/j e half-space DtN operator defined by 



(4.1), 



G £(i7" 1 / 2 (E H ),ff- 1 / 2 (E )) is the restriction operator defined by fi~4) 



E 27T / 3 G £(if 2 \ 3 (£o), .ff^Yg^E')) is the extension operator by symmetry 
defined in Section \3. S\ 



1/2/ 

Proof. Let <^> G ^/^^(E 1 ). From the definition of D 2lT / 3 , the functions w e (</>) |n« 



and u H (D 27r /3<fi) satisfy the half-space Helmholtz problem (4.2 ) with the same Dirich- 
let condition on H H , namely ip — D 27T / 3 <p. The uniqueness of the solution of this 
problem implies that u e (cf>) \qh = u (D 27r / 3 4>) and in particular, the traces of their 
normal derivatives on E° coincide, yielding 



^2ir/3 



du e ((t>) 



(R H o A H o D 2n/3 ) 



where R 2v / 3 is the restriction operator by symmetry defined in Section 3.2 This 
relation proves the first part of the theorem. For the second part, we just use that 
E 27r / 3 is the inverse of R 27r / 3 - D 

In section [3J we explain how to compute the half-space DtN operator with the help of 
an adapted version of the Floquct-Bloch transform defined in|5.1.1|and the resolution 



of a family of half-space problems (4.2 1 with k— quasiperiodic boundary conditions. 



Section[6]deals with the characterization of the DtD operator D 27T / 3 . The computation 



of D 2v / 3 a priori requires to compute the solutions u e of the exterior problem (2.2 1 
defined in an unbounded domain. We explain, using the half-space problem and the 
properties of the problem, how to obtain a characterization of this operator which 
avoid the solution of the exterior problem. 

5. Characterization of the half-space DtN operator In this section, we 
tackle the Dirichlet half-space problem. In other words, for any <\> G i? 1 / 2 (E if ) we 
want to compute the solution u H (4>) in H 1 (A, tt H ) of 



Au H ((f)) + pu H ((/)) =0, infi H , 
u H (cj)) = (/>, on E H , 
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(see Figure 2.1 for notations). We will deduce a characterization of the half-space 



DtN operator A H 

du H (tj>) 



A' 



where v H is the exterior normal to il H . 

Remark 5.1. The half-space problem is not only interesting as a step of our ap- 
proach to analyze transparent boundary conditions for locally perturbed hexagonal pe- 
riodic media. Indeed, it also appears naturally in transmission problems between a 
homogeneous medium and a hexagonal periodic one (see where such problems are 
considered for the case of one dimensional type periodic media and J18\[ for the case 
of square lattices.) 



We develop a method for computing the solution of (4.2) and the operator A H , by 
adapting the method developed in [T71HH]. In these works, the half-space solution and 
then the half-space DtN operator are computed using the Floquet Bloch transform. 
More precisely, this is done via the resolution of a family of waveguide problems 
with quasiperiodic conditions, each waveguide solution being computed thanks to the 
resolution of elementary cell problems and a stationary Riccati equation. In our case, 
this approach cannot be directly transposed since 

• the waveguide boundary would not correspond to a physical boundary (see 



the shaded domain of Figure 5.1 ); 
• but most importantly, computing D^fi would be much more intricate (see 



Remark 6.4 1 



Instead, the half-space problem is handled by solving a family of half-space problems 
with k— quasiperiodic boundary conditions. 

5.1. The half-space problem: from arbitrary data to quasiperiodic data 

5.1.1. The Floquet-Bloch (FB) transform and its properties Following 
[T7l [28] , we recall below the definition of the FB transform and state without proof 
its main properties. 

Set K = Ix (-tt/L,tt/L). The FB transform of period L (here L = y/3d) is de- 
fined by 

T: C§°(R) L 2 QP (K) 

where Lq P (K) is the set of functions / G Lf oc such that for any k E (—ir/L,ir/L), 

f(-,k) is fc— quasi-periodic, that means f(y + nL,k) = f(y,k)e* . We equip this 
space by the norm of L 2 (Ko) where Ko = (— L/2,L/2) x (—ir/L,ir/L). 

The operator T can be extended as an isometry between L 2 (K) and Lq P (K): 

V0, i> G i 2 (M), (T<f>, ^) L 2 Qp(K) = (4>,ib) L 2 (M) . 

This transform is a privileged tool for the study of PDE with periodic coefficients 
because it commutes with 
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• any differential operator ; 

• the multiplication by any periodic function with period L. 
The FB transform is invertible and the inversion formula is given by 



<t>(y) = 



tt/L 



J-(j)(y; k)dk. 



(5.1) 



-tt/L 



Next we define the partial FB transform in the y-direction in the half-space £l H : 



L 2 (n H ) - 

u(x,y) i-> T y u(x,y;k) 



.ZvQp ( X 



7T 7T 

V L 



with 



Vz, {F y u){x, •;•) = F[u(x, ■)] 



and Lq P (Q h x (— n/L, tt/X)) is the set of square integrable functions /, locally in 
the y— direction, such that for any k G (—n/L,ir/L), /(•, k) is k— quasi-periodic in 
the y— direction, that means f(-,y + nL,k) = f(-,y,k)e mkL . We equip this space 
by the norm of L 2 (ft w x (-ir/L, vr/L)) where Q w — fl H n{y€ (-i/2,£/2)} (see 
Figure 5.1 ). 




Fig. 5.1. The domain Q w and its left boundary EL- 



It is easy to see that the partial FB transform J- y defines an isomorphism from L 2 (fl H ) 
into L 2 QP (fl H x (-tt/L,tt/L)). 

Now we want to know how the Floquet Bloch transform is defined or can be ex- 
tended to the other functional spaces appearing in our study (H 1 ^ 11 , A), i? 1 / 2 (S i/ ) 
and iJ _1 / 2 (E i/ )). To make a rigorous presentation, we need to introduce, k being a 
parameter between —n/L and iv/L, the so called fc-quasiperiodic extension operator 

E^ P G C (L 2 (n w ),L 2 untf (n H )) defined by 

Vu G L 2 (Q W ), Vq G Z, V(x, y) G Q W , E® P u{x, y + qL) = u(x, y) e lqkL , 
where L 2 ini ^(Q, H ) is the normed space defined by 

L 2 vnif (n H ) := \u G L 2 oc (^), sup / M 2 < +ool . 
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A natural functional space which appears is the set of locally L 2 , k— quasiperiodic 
functions defined in tt H , denoted L 2 (il H ) and characterized by 

We can introduce the corresponding fc-quasiperiodic restriction operator 

RQ p ec(Ll(n H ),L 2 (n w )) 

defined by 

Vu fe e Ll(n H ), R^ P u k := u k \ nW . 

Noting that 

L'{n") = B$ p {Ltta*)) t 

it is easy to see that L 2 (Q H ) and L 2 (fl w ) are isomorphic. We can then consider 
L 2 (il H ) as a Hilbert space when it is endowed with the inner product of L 2 {Q W ). 
More precisely, we define the scalar product on L 2 (il H ) as follows 

v{u k ,v k ) e L 2 (n H ) 2 , (u k ,v k ) L 2 (nH) = (R% p u k ,R% p v k ) w . 

Next we define smooth quasiperiodic functions in n, H 

C^(n H ) = {ueC°°(n H ), u(x,y + L)=u{x,y)e lkL : \/{x,y) G tt H } . 
and smooth quasiperiodic functions in Sl w 

cr(n w ) = R% p (cr(n ff )) 

Let Hl(fl w ) (resp. Hl(A,Ct w )) be the closure of C^{n w ) in H l {Vt w ) (resp. in 
ff^A, fl w )) equipped with the norm of H 1 ^) (resp. the norm of if 1 (A, Cl w )) and 
let Hl(n H ) (resp. i^(A,ft H )) be defined by 

Ht(n») = E% p (Hl(n w )) ( resp. /T*(A,n*) = (^(A,0^))) , 

that we also equip with the norm of H 1 (Sl w ) (resp. the norm of H 1 (A,Q W )). 
Remark 5.2. The functions of Hl (n H ) are nothing but the k— quasiperiodic exten- 
sions of functions in Hl(il w ). Therefore they are H 1 in any horizontal strip but not 
in the vertical ones. 
Let us denote (see Figure pTT] ) 

14 :=z H ntW y= s°) . 

and let us define H^ 2 CE,Q a ) by 

equipped with the graph norm. We define the space H^ 2 (E H ) by k— quasiperiodic 
extension of functions of h1^ 2 (J£q ): 
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that we equip with the norm of LT/J 2 : (Sg ). 

The space /f^ 2 (£oo) is a dense subspace of H 1 ^ 2 (T,l ) and the embedding from 
i?y 2 (Sq ) onto H 1 ^ 2 (Y.q Q ) is continuous. We can then dehne the dual space of 

H^ 2 (Eq ) that we call H^ 1 ' 2 ^^). According to Green's formula, we can show 
that 

du 
dv 



H -l/2 (y i s 



,«eff t 1 (A,!l H ') 



To define the space of /c-extension of functions in H, 1 ^ 2 (Sg ), we need to define the 
extension operator E? P on H// 1 ^ 2 (Sqo) m a wea k sense. Actually, this can be done 
by duality by setting for all (ip Q ,<t> k ) € # fc ~ 1/2 (^ ) x ^ 1/2 (S H ) 

Finally, we define iJ- 1/2 (£ H ) - £« p (^ 1/2 (E^ )) . 

We can now state the following results. 

Theorem 5.1. J- y is an isomorphism from X H onto 

X QP := {u E L 2 (-tt/L, tt/L; X H ) \ for a. e. k E (-tt/L, tt/L), u(-:k) E Xf } , 

j-Tv/L 

equipped with the norm \\u\\ 2 x = / fc)!) 2 ^ <ifc where 

J-n/L k 

• A H =JJ 1 (n ff ), X QP = H l QP (Q, H x (-tt/L, tt/L)) and X^ = Hl(fl H ); 

• X H =H 1 (A,9, H ), A QP = i^ P (A;f7 ff x (-tt/L, tt/L)) and 
XX = Hl(A,n H ); 

. A ff = i/ 1 /2( S ff) ; Xgp = #V p 2 ( E H x {-tc/L^IL)) and X« = # fc 1/2 (£ ff ). 
Finally, we can extend by duality the definition of T y to the space H~ 1 ^ 2 (Y, H ) intro- 
ducing the dual of H^ 2 (Z H x (-tt/L, tt/L)) 

H Q l' 2 (H H x(-tt/L,tt/L)) := 

G L 2 (-vr/L, tt/L; A) | Vfc G (-tt/L, tt/L), fc) G ff,7 1/2 (£ H ) } , 

whereA = £ QP (7i4 /2 (l4)" 



According to Theorem |5.1[ the partial FB transform JF y defines an isomorphism from 
H 1 / 2 (Y, H ) onto (Y* h x (— tt/L, tt/L)). Using the Riesz representation theorem, 
the FB transform can then be extended by duality as an isomorphism from 77 _1 / 2 (S fr ) 
onto H Q p /2 (T, H x (-tt/L, tt/L)) see [in] for more details. 

5.1.2. Application to the half-space problem The above results imply 
that for almost every k in (—tt/L, tt/L), we have 

G H 1 / 2 ^), fa := k) e Hl /2 (^ H ) 

W G H- 1 / 2 ^ 11 ), i> k :=F y ^(-;k) E H- 1/2 (Z H ) 

Vu H E ff^A.ft*), u% :=T y u H (-,k) E Hl(A,n H ) 
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The following theorem is a direct consequence of the properties of the FB transform 
given in the previous section. 

Theorem 5.2. Let u H (cj)) be the solution of problem (V H ) (see equation ( |4.2[ ) ). 
For every k £ (—n/L,n/L), u^(</>fc) :— J- y iu H {(j))^ (•; k) is the unique solution 
Hl(A,fl H ) of the half-space problem with the k-quasiper iodic boundary condition 
J- y (j)(-;k). Using the inversion formula (5.1|, we can recover u H (4>) in the 



whole domain Q H for any Dirichlet condition <f> £W-^ 2 (T, H ): 



= / u^{<j> k )dk. (5.2) 
V 2tt J_^ /l 

Then we can show the following theorem which expresses that the half-space DtN 

operator A H can be described in terms of a family of "quasiperiodic" half-space DtN 
operators. 

Theorem 5.3. The half-space DtN operator A H is given by: 

e HV 2 (E H ), A H ct>=J^r L Af0 k dk, (5.3) 
V 2tt J_ w/l 



where A^ is k— quasiperiodic half-space DtN operator, defined by 

h 2 d ^kGk) 



A 



(5.4) 



According to relations (5.2) (resp. (5.3)), the solution of the half-space problem 
(resp. the half-space DtN operator) for arbitrary boundary data (f> is obtained by 
superposing the corresponding solutions (resp. DtN operator) for quasiperiodic data. 
The next subsection is devoted to solving such problems. 

5.2. Solution of the half-space problem for quasiperiodic boundary data 

Let k be in (— n/L, tt/L), we explain here how to compute the solution of (V ) (see 

=]> k e Ht /2 (E H ). We have seen in the 
and ! 



(4.2)) for k~ quasiperiodic boundary data (f> 
previous sections (see in particular Theorems 



5.2 



5.31 that for any 4> k £ H\' (T, H ) 



fO) admits a unique solution u%Q> k ) € Hl(A,tt H ) and Af 4> k E H k 1 (£ H ) 



This half-space problem is in some sense the counterpart of the waveguide prob- 
lem with fc— quasiperiodic conditions used in [TTHH] to determine the half-space DtN 
operator for the case of a square periodicity cell. In jTTl QjJ] , the basic tools are the 
resolution of elementary cell problems and a stationary Riccati equation whose solu- 
tion is a so-called propagation operator. 

Let us begin with some notation. Let Coo be a periodicity cell whose boundary 
meets the vertical boundary T, H . Given p E N, q E Z, we introduce the vector 
V pq = pei + qe2 (See Assumption (Al) in Section ljand Figure 1.2 for the definition 



of the directions of periodicity ei and e 2 ). The eel 
then be defined by translation of the reference cell (see Figure 2.3 ) 



C pq of the periodic half-space can 



Vp E N, q E 



Cpq — CqO 



V 



PQ- 



We will denote by J7 p the vertical "strip" containing the cell C p o: 

Qp Cpq . 
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In the following, for a cell of periodicity C pq , we introduce the oriented boundaries 
described in Figure [572) 




Fig. 5.2. Notations used for a periodicity cell 



Finally, using the fc-quasiperiodic restriction operator Rj: and the fc-quasiperiodic 
extension operator E® P defined in Section 5.1.1 we recall that H]J 2 {Y> H ) is isomor- 



phic to Hl /2 (Z l 00 ). 

5.2.1. The propagation operator We can now introduce the operator Vk 
defined by 



where u^(E^ P <f>1) is the unique solution of (V H ) (see (4.2)) with boundary condition 
the k— quasiperiodic extension of c/f>°: E^ P (fi®. 



< 



Fig. 5.3. Description of the propagation operator 



Remark 5.3. For the sake of simplicity, we will often identify throughout the paper 
functional spaces of functions acting on with the same spaces acting on Sqo- 

Typically, iZ^ 2 (E' ) for arbitrary (p, q) G N x Z will be identified with H^ 2 (T,q ). 

The operator Vk, considered now as a bounded linear operator from H^ 2 (Eq ) onto 
itself (see Remark 5.3 1, is called a "propagation operator". Indeed, the solution can 
be reconstructed in any cell of the medium from its values in the reference cell Coo 
and the knowledge of Vk, as shown in the next result. 

Theorem 5.4. For any k— quasiperiodic condition <f>k £ H^ 2 (T, H ), the solution 
^k&k) of p~2| is given by 

VpeN, q e Z, 



u H k { 



e^u^ErTiRT, 



(5.5) 
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Proof. First of all, it is clear that the periodicity and the well-posedness of the 
problem in the vertical direction implies that for any q € Z, we have 



Coo 



(5.6) 



In the horizontal direction, due to the periodicity and the well-posedness of problem 
(4.2), we note that the value of u on the cell C10 is related to the one on the cell Cqo 



as follows 



(5.7) 



By induction, one easily gets that for all p € N there holds 



QP ■ 



(5.8) 



Combining relations (5.6) and (5.8) yields the claimed result (5.5). □ 



The next result collects some useful properties of the propagation operator Vk- 
Corollary 5.5. The operator V k € £(i/fe /2 (Sg )) is a compact operator with spectral 
radius strictly less than one. 

Proof. The compactness of V k € C(H^ 2 (T,q )) follows easily from interior regu- 
larity and Sobolev compactness embedding arguments. 



By definition of the space Hl(A, £l H ) (see Remark 5.2), (<fi k ) satisfies in particular 



for any (f> k G # fc 1/2 (£ H ) 



/ |£f(0 fe )| 2 <+«). 



Due to the fc-quasiperiodicity of (cf> k ), the above relation is equivalent to 



\u" (?fe)| 2 < +™ 



Moreover, using Theorem 5.4 we have 



%$k)\ 2 = E / is* (&)i s 



E/ \^(e, 

tm « Cnn 



Consequently, if A is eigenvalue of V k and ip an associated eigenvector, combining the 
last two relations for 4> k = ip shows that 




|nf M| 2 < +00 



from which we get that 



IAI <1. 
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□ 

According to Theorem |5.4| the solution of the half-space problem (4.2) for a k— quasi- 
periodic Dirichlet condition is completely determined on the whole domain fl H as 
soon as it is known on the reference periodicity cell Coo, provided the propagation 
operator Vk is also known. 

5.2.2. Elementary problems 

Like in [TTJ [TH] , introducing elementary problems allows us to restrict the half- 
space problem for quasiperiodic boundary data to the determination of the propaga- 
tion operator Vk- 

More precisely, given 0° € H^ 2 (12q Q ), let us introduce the solutions of the follow- 
ing problems set in the vertical strip Q = U q ez ^-o<? 



E^^k) "= -^fc(^>^o) is the unique solution of (see Figure 5.4 1 




0, in fi , 

on Ego, (5.9) 



on Efo, 



• Ek(</>1) £ HI(A,£Iq) is the unique solution of (see Figure 5.4) 

A££($) = 0, infio, 

E r M)=0, "nS;, (5.10) 

E r k (<t>l)=<t>l ™Ef . 

Let us emphasize that the boundary conditions are described only for the two left 
lateral sides of Coo and Cio, as the boundary condition on Eq ? (resp. on S| ? ) for 
q € Z* follows directly from the one on Sq (resp. on E^ ) by k— quasiperiodicity. 




Fig. 5.4. The two strip problems for Ef (on the left) and ET (on the right) 



We also introduce the elementary cell solutions 

4($) = EM)\ Cqo and eM) = ^(/ fc )| Coo • 

Conversely, it is clear that £7£(</>°) and ££(<^) are uniquely determined by the above 
elementary cell solutions due to their quasiperiodicity. Therefore, in practice, one 
simply needs to solve the corresponding cell problems set in the reference periodicity 
cell Cqo- One might thus wonder why we have introduced the strip problems and 
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not directly the cell problems. In fact, it turns out that viewing the elementary cell 
solutions as restrictions of the strip problems leads to more compact and simpler ex- 
pressions in the following. 

The main advantage of these elementary problems lies in the fact that, by linear- 



ity of (4.2 ), one has 



and then in the reference cell 



4($) + e*(P*$). 



(5.11) 



5.2.3. The Ricatti equation for the determination of the propagation 
operator 

Assuming the elementary cell solutions are known, it remains to determine the 
propagation operator V k . To this end, we use the same strategy as in [T7]. In short, 
the equation characterizing the propagation operator Vk is obtained by writing the 
continuity of the normal derivative of uff across each boundary which is ensured 
by Theorem |5.4| 

To this end, we introduce four local DtN operators associated to the elementary 



problems (5.9) and (5.10|. We refer the reader to Section 5.1.1 for the definition of 
the spaces Hj 2 (Ys l m ) and H k 1 ^ 2 (Sq ). 

Definition 5.6. We introduce the following local DtN operators 

V h i G C (i^ 1/2 (s4),i^ 1/2 (s4)) , i,j G {i,r}, 

where for all <jP k G H^ 2 (Eq ) : 



7. Kfk) ' v 



I, 1 7. - x ^I ' Vk) ' v 



where v is the outgoing unit normal to C o . 




Fig. 5.5. The four local DtN operators 



The characterization of the operator Vk is then given by the following result. 
Theorem 5.7. The propagation operator Vk is the unique compact operator of 
£{Hl (Sqo)) with spectral radius strictly less than 1 solution of the stationary Riccati 
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equation 



r k re r 2 k + (T k " + n r )v k + n 



-ff 



(5.12) 



Proof. The proof involves two steps: we first show that V k satisfies (5.12) and 



then we prove that (5.12| has a unique solution with spectral radius strictly less than 1. 



Step 1: V k satisfies (5.12 1. First of all, according to Corollary 5.5 V k is an oper- 
ator in C{h]/ 2 (Sqo)) w hh spectral radius strictly less than 1. Moreover, using the 
continuity of the normal derivative of u k , we obtain using equations (5.7) and (5.11 ), 



that for any <f>° k £ H L k /z (T, e 0Q ) 



which gives 



•7TW2) = -T^iv^D - r k r \vU\). 



Since $ £ iJ fc 1/2 (^ 



is arbitrary, we get 
r k er + (T k rr + T k u )V k +T k re V 2 k 



0. 



which is exactly (5.12) 



Step 2: Uniqueness for (5.12). This result follows immediately from the unique- 



ness of the solution of (4.2). Indeed, assume that V £ £(H^ 2 (Y,q )) has a spectral 



radius strictly less than 1 and satisfies 



rrv 2 + (T k 



T k rr )v + rt = o. 



1/2 



(T, H ) , the function v defined in 



Then, one can easily show that for any <f> k £ H k 
each cell C pq , pGN, q £ Z, by 

v \c pq = e"? fe M(^ Rg P ?*) + e r k (V p+1 R2 P 

belongs to Hl(A,fl H ) (thanks to the definition of e e k and e k and using the fact that 
V is a solution of the Riccati equation with spectral radius strictly less than 1) and 



solves (4.2). Then, the uniqueness of the solution of (4.2 1 in Hl(A,H, H ) implies that 

v = u k {4> k ), 

and by definition of v and V k , we deduce that V = V k . □ 



By Theorem 5.4 and expression (5.11), solving the elementary problems ( 5.9 (-( 5.10 1 



and the Riccati equation ( 5.12 ) allows us to reconstruct cell by cell the unique solution 



of (V H ) (see (4.2)) in the case of quasiperiodic boundary data. 



Finally, we deduce from the above analysis the expression of the DtN operator for 
quasiperiodic boundary condition: 

Proposition 5.8 (DtN operator for quasiperiodic boundary data). Suppose that 
the data <f> k £ ii\J 2 (T, H ) , the DtN operator is given by 



E 



Qp r k "Rr- 



E 



Qp r k re r k Rf p , 



(5.13) 



Transparent boundary conditions for locally perturbed infinite hexagonal periodic media 23 



5.2.4. Additional tools We conclude this subsection by introducing four local 
Dirichlet-to-Dirichlet (DtD) operators associated to the cell Coo that will be needed 
in the sequel. 

Definition 5.9. We define the local DtD operators by setting for all <jP k e H^ 2 (Eq ): 



EM) 



and 
and 



v r k ~<i>l 



mi) 

E r M) 



Identifying Sq , r oo and Fq Q , the above operators will be considered as bounded linear 
operators from iJ^ 2 (Soo) onto i? 1//2 (E u0 ). 







Fig. 5.6. The four local DtD operators 



Using Theorem 5.4 and expression (5.11) we obtain the following result 



Corollary 5.10. For any 4> k e HjJ^^E, 11 ), for any pgN and q G Z, we have 



— e mkL 



V[ + VtR Q k P ^ + V r k + Vi +1 R Q k P > 



u H k { 



iqkL 



Vl-ViR Q k P ^ 



6. Characterization of the DtD operator This section is devoted to the 
determination of the DtD operator D 27T / 3 ■ We first show in Section 16.11 that it solves 
an affine equation involving an operator D H associated with the half-space problem. 
Using the FB transform, we derive a semi-analytic expression for D H in Section 
6.2| Finally, we deduce in Section |6.3| an equivalent integral formulation of the affine 
equation which is more suitable for future numerical approximation. 

6.1. The affine equation 



First of all let us introduce some useful notation (see Figure 6.1 ): S ff = S US°UI] + , 
^2tt/3 = ®27r/3S~, ^2?r/3 = ^-2tt/3^ + ; where the center of the rotations 0±2tt/3 is 
the center of f2* . Let us remind that £l H is the half-space at the right of Yj H . 



Let us recall that D 27r /3 is defined by (see Figure 2.2) 

D 27r/3 : JJV2(E<) _> ¥/ 2 (^) 

<j> I > U e ((j>) | S H . 

We first remark that D 27T / 3 belongs to the affine space: 
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Fig. 6.1. Symmetry axis 



Let us introduce, D H , called the half-space DtD operator, defined by 







s _ = u h (iJj) 






s o = " H (V0| 






s+ = u H {^)\ 



S 2i/3 



where we have identified S and E 2ir ^ 3 , S + and S^. , 3 taking into account the direc 



tions shown in Figure 6.1 



Remark 6.1. The range of D H is included in ff 1 / 2 (S ff ) because for any ip € 
iJ 1/,2 (S ff ), D H ij} is nothing but the trace of the H 1 function u H (ip) on the broken 
line ^2tt/3 U S° U E^. , 3 (identifying Y> H with S^. / 3 U S° U £J, 3y ). 
We have then the following fundamental theorem 

Theorem 6.1. TTie operator D 2 ^/ 3 is the unique solution of the problem: 

Find D G £ S o suc/i i/ioi D = D H oD. (6.1) 



Remark 6.2. ./Vote i/ia£ since £s° is an affine space, the problem (6.1) is an afjine 
problem, even though the equation is linear. 

Proof. Existence: We prove that the operator D 27r /3 is a solution of (6.1). We 
have already seen that 



(6.2) 



Moreover, since <j) is in H^ 2 , 3 CE l ), u e ((f>) is in .ffj , 3 (fi e ). In particular, 
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^2ir/3 = ®27r/3^ which implies u e (0)| s - = u e (</>)|s- and then 



i H (D 2 „ /3 <l>)\ s - =u e (0)| s - using @, 

2tv /3 



S + = 62^/3!]^.^ which implies u e (0)| s + = u e (</>)| E + and then u e (<j>)\-£+ = 



M ff (-D27r/30)I E + using ( [672 1 



2ir/3 



Using the definition of D H , we obtain that D 27T /3 solves (6.1 1. 



Uniqueness: Let D be an operator from H^ 3 (Y^ 1 ) into i/ 1 / 2 (S H ) such that for all 

1/2 

4> G iJ 2 ^y 3 (S l ), -D^lso = and which satisfies 



D H o D - D = 0. 

12 

r/3 



(6.3) 



We prove that D = 0. Let <p e H^ /3 (T^) and w = u H (D<j>) defined in We have 
in particular 



v \xo =u H {D4>)\vo =Dcf>\ & = 0. 

Now we build a function in the half-space £l 27T / 3 — < ^ > 2w/3^ H by: v 2w /3 — vq{Q 27T / 3 x). 
By a classical argument, since Vq is solution of (4.2) in it is clear that ^2^/3 is 
solution of 



Aw 2 7r/3 + P V 2*/3 = °: in ^2^/3; 



while vqIY, = implies ^2^/3 e 2 , r / ! £ = 0. We are going to show that Vq and ^2^/3 
coincide in fl H D (Q 27T / 3 il H ). The difference d 2v / 3 — v o ~ u 2tt/3 satisfies 



Ad 27l/3 + pd 27T/3 = o 7 in r> ff n (e 27r/3 n H ), 

with boundary condition: 

^2^/3 Ie+ = u |s+ - "27T/3 ls+ = V Q ls+ ~ v 



(6.4) 



^2ir/3 



^21/3 



£- - ^2^/3 

2ir/3 ' 



2tt/3 



^0 |e- 



since £+ = e2 W / 3 S+. /3 , S 2Tr/3 = 6 23r / 3 £ and t; 27r/3 = v (e 2 ^ /3 ). Using the defini- 
tion of vq , we have 



<Ws|s+ = « H (/ty)| E + -U H (D(/)) 



^2^/3 



which gives using the definitions of u ff and D H 



and then by (6.3) 



<W 3 |e+ =^|e+ -£ H °^| s+ , 



^27r/3 1 2+ — 0. 



(6.5) 



2G 



Besse, Coatleven, Fliss, Lacroix- Violet, Ramdani 



In the same way, we have 



2tt/3 



u H {D(f>) 



2ir/3 



u H {Dcp) 



and then, using again the definition of u H and D H 



2 2tt/3 



Dcf) 



D H oD<t> 



0. 



(6.6) 



A uniqueness argument for the Dirichlet problem ( 6.4 ), (6.5), (6.6) yields ^2^/3 = v o 

in n H n (e 27r/3 n ff ). 

With the same argument, we can construct three solutions of the Helmholtz equation 
v 0i v 2tt/3i v 4tt/3 respectively in the domains Q H , Q 2 -k/3^ H an d @4 7r /3^ i? and which 
coincide in the domains where they are jointly defined: 

• v = v 2rr/3 in n H n (e 27r/3 ft H ), 

• W27T/3 = V4tt/3 in (02tt/3^ H ) R ( Q ^/3^ H ), 

• v i7r/3 = v in (e 47r / 3 ^ H ) n n H . 

Thus we can construct a function C/ e e iJ 1 (fi e ) defined in J7 e by 

^ |e w ,n« = VfeTr/3, * e {0, 2, 4}, 

so that C/ e is an if 1 -function which satisfies 

At/ e + p[/ e = 0, inft e , 

with homogeneous Dirichlet condition on £*. From the uniqueness for the exterior 
problem we get U e = in fi e and therefore Dtp\^i = f7 e |s« = 0. This concludes the 
proof of Theorem 6.1 □ 

6.2. Characterization of the half-space DtD operator 

Let us recall the definition of D H : 

D H : H X / 2 (Y, H ) — > H 1 / 2 ^") 





E- 


= U H (V>) 


2; 


D h i[> 


E° 




£° 


D H ip 


S + 


EE U H (lp) 


2; 



277/3 



where we have identified £ and S 97r ^ 3 , S + and S^. , 3 taking into account the direc- 
tions shown in Figure [6~Tj 



Thanks to the results of Section [5] we can give a semi-analytic expression for the 
DtD operator D H . 

Proposition 6.2. For any ip in if x /2(s ff ) and for any £ in (—tt/L, tt/L), we have 



F y {D H <P)(;0 = ~ e <L (?{+ + V?V k ) (I V k e* L ) 1 B$ p $ k dk 

llT J -tt/L 



L 

2k 

L 
2^ 

L 
2^ 



-n/ 

-tt/L 
tt/L 

1 

-tt/L 



i>k |eo dk 



(6.7) 



-»(*!+£)£ 



{V{-+Vl-V k ) (l-V.e-^+^y 1 R^ k 



dk, 
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where tpf. = J^yipl^k), the local DtD operators and VI* are given in Definition 
\5.9\ and Vk is the propagation operator defined in Section \5.2.1\ 



Remark 6.3. Let us remark that the subscript k of the operators in formula (6.7 1 
means that they actually depend on the variable k. 

Proof. Let k £ (— n/L, n/L), we express first D H %pk for k— quasiperiodic bound- 
ary data ipk £ (E H ). By definition of D H , we are then interested in expressing 



remark that 



Uk(ipk) |s° and u"{4>k) 



Using notation of Section 



5.2 



let us 



S 2tt/3 — U ^nO' S 2*73 ~~ U ^n,-(n+l) > 



n=-\-oo n=0 

taking into account the directions shown in Figures |6.1| and |5.2| 



Using Corollary |5.f0[ we have 
• on Yj~ for all n £ N, 

ZTT / 3 7 7 



and D H i[>k 
on S° 



Ukfyk) p+ wnere we have denoted £„ = 0_ 2 7r/3r^o; 



-D^^fe |e° = «f Ws) Is" = |e° ; 



• on DJiyg, for all neN, 



where £+ = 5 _/ s r 



27T/3 1 n> _(„ + i) • 



For any data ^ 6 i/ 1 / 2 (E ff ), we use the FB transform and formula (5.2) to obtain 
the expression "piece by piece" of D H ip from -ipk = Fyipi', k) 
• on XT = U„ e NS^, for all n £ N, 



-Tf/L 



on £ u 



r/L 



Ie° dfc ; 



-7V/L 



on E + = U neN £„ , for all n £ N, 



7r/L 



7T/L 



28 



Besse, Coatleven, Fliss, Lacroix- Violet, Ramdani 



We finally apply the FB transform in the y— direction to D H ip: 

o 



n=0 



^+ e 



-i£(n+l)L 



By inverting the integral over (—tt/L, tt/L) and the sum over n, we are led to use the 
following formula: 



n<EN 



n e ±!(n+l)Ci 



(6.8) 



with ( — fc or £ = fc + £. Note that this is possible as for every k, Vk is compact 
with spectral radius strictly less than 1. Actually, we could prove like in [T7] that for 
Pb > (defined in (1.5)), the spectral radius p(Vk) of Vk is uniformly bounded in k, 
by a constant that is strictly less than 1: 



3C> 0, Vfc G |) , p(7> fc )<e- 



Cp 6 



The property 



lim H^ll 1 /" = p(p fc ) 

n— f +oo 



for the norm of £(L 2 (Eq )) ([il]) implies that for some a €]e Cpi> , 1[, j large enough 
we have for all fc 



which yields the absolute convergence of the series (6.8). Therefore, for each £, X — 
T'ke ±l ^ L is invertible and the sum (6.8 ) converges uniformly in the norm of £(L 2 (Eq )). 



Exchanging the order of the integral and the sum is then possible. □ 
Remark 6.4. Our choice of solving a family of k— quasiperiodic half-space problems 
instead of a family of k— quasiperiodic half waveguide problems simplifies the expres- 
sion of D H . Indeed, expressing u H (<fi) on S^, 3 is much easier by doing so. 

6.3. Towards the integral equation 

Let us give now more precisions about the resolution of the affine equation (6.1 ). 



Since the operator D H is characterized via its FB transform, it makes sense to refor 



mulate (6.1) using the FB transform. 



1/2 

Corollary 6.3. For any cj> £ / 3 (S l ), the function 



is the unique solution to the following problem 



Find $ E H^p (T, H 



V L 



such that 



(u) 



■k/L 



f ' K H (k,£)^{-,k)dk = 0, 

J-tt/L 



(6.9) 



1p(; k) | E o dk = (j) | s o, 



-tt/L 
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where K H (k,£) G C (ff^E^), iI £ 1/2 (I! H )) is given by 



K H (k, = — e< L (V{+ + V r k + V k ) (1 - V k e 



, R QP 



L 

2^ 
L 

— < 
2tt 



R 



H 



<L\ 1 nQP 



with R H the restriction operator from Yi H on E° defined in Section [^J Relation 
( 6.9 is the FB transform of Equation (6.1| while relation (6.9)-(ii) expresses in 
terms of the FB- variables the condition: 



D 



2tt/3<P |S° 



From a practical point of view, it seems to us that it is easier to solve (6.9) instead 



of (6.1). The advantage is to replace the discretization of an infinite set by the 



discretization of a compact set. 

7. Summary and conclusion In this paper, we proposed a method to solve 
scattering problems in infinite hexagonal periodic media containing local perturba- 
tions. We have computed the Dirichlet-to-Neumann map on the boundary of a 
bounded domain which respects the hexagonal geometry. By doing so, the initial 
problem is reduced to the solution of a boundray value problem set in a cell con- 
taining the defect. The computation of the DtN map A is based on a factoriz ation 
through a half-space DtN operator A H and a DtD operator D 2 ^/^ (see Theorem 4.1 ). 
We sum up below the main steps to be followed in order to compute A</>, for a given 
boundary data with hexagonal symmetry 
1. Pre-processing steps 



For all k G (-ir/L,ir/L) (a) Solve the cell problems ( 5.9 1-( 5.10) 



(a 
(b) 



Compute the local DtN operators 



(see Definition (5.6)) 



(c) Compute the local DtD operators 



(see Definition (5.9)) 



(d) Solve the stationary Riccati equation (5.12) 



(e) Compute A^ thanks to relation (5.13) 



Solve integral equation (6.9) using the steps (l)-(c) and (l)-(d) to obtain 
F v (p. 



1> 



27r/3 



2tt/3V 



3. Using relation (5.3), 



A< 



7T/L 



-7t/£ 



AT^2V3('> fc ) dk 



The pre-processing steps only involve the solution of elementary problems set on 
bounded domains, for each values of k. These steps can be easily parallelized as the 



problems for different values of k are decoupled. Solving the integral equation (6.9) 



constitutes the main difficulty in implementing the algorithm. The discretization and 
numerical investigation of the problem will be considered in a forthcoming article. 
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